J /if} suppression in Pb+Pb collisions and pr broadening 
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We have analysed the NA50 data, on the centrality depen- 
dence of pt broadening of J/tp's, in Pb+Pb collisions, at the 
CERN-SPS. The data were analysed in a QCD based model, 
where J/ip's are suppressed in 'nuclear' medium. Without 
any free parameter, the model could explain the NA50 pr 
broadening data. The data were also analysed in a QGP 
based threshold model, where J/ip suppression is 100% above 
a critical density. The QGP based model could not explain 
the NA50 pr broadening data. We have also predicted the 
centrality dependence of J/ip suppression and pr broadening 
at RHIC energy. Both the models, the QGP based threshold 
model and the QCD based nuclear absorption model, predict 
Pt broadening very close to each other. 



Since the prediction by Matsui and Satz [1] that bind- 
ing of a cc pair into a J /if) meson will be hindered in 
quark-gluon plasma (QGP), J/ip suppression is recog- 
nized as an important tool for the identification of the 
possible phase transition from confined to deconfmcd 
matter. NA50 collaboration measured centrality depen- 
dence of J/ip suppression in 158 A GeV Pb+Pb collisions 
[2]. They observed suppression well beyond the stan- 
dard nuclear absorption model. Initially the data were 
interpreted in terms of successive melting of charmonium 
states in QGP [2]. However, later, it was realized that 
the data could be explained in a variety of models with 
or without QGP [3-7]. What it more intriguing is that 
the predicted centrality dependence of J/ ip at RHIC en- 
ergy in a model without QGP, matches with the model 
prediction with QGP [6]. It appears that the J/ip sup- 
pression may not be a good signal for the deconfining 
phase transition. 

Apart from the centrality dependence of J/ip sup- 
pression, NA50 collaboration also presented data on the 
transverse energy (centrality) dependence of pt broad- 
ening of J/f/Vs [8], which did not receive much attention. 
Kharzeev et al [9] suggested that px broadening of J/ip 
can distinguish between QGP and nuclear matter. They 
argued that in a nuclear matter, px broadening will sat- 
urate at large Et- In contrast, in a QGP, pr broadening 
will visibly decrease. They used conventional Glauber 
model of nuclear absorption, which we know, could not 
explain the NA50 data on the centrality dependence of 
J/ip over Drell-Yan ratio. The centrality dependence of 
J / ip suppression is well explained in the 'unconventional' 
QCD based nuclear absorption model [6]. In the present 
letter we have analysed the NA50 pt broadening data in 



the model. Without any free parameter, the model could 
explain the data. For comparison purpose, we have anal- 
ysed the data also in the QGP based threshold model [3] . 
The model fails to fit the data. We have also predicted 
Pt broadening at RHIC energy. Interestingly, both the 
models predict very similar pt broadening at RHIC. It 
seems that, like the J/ip suppression, centrality depen- 
dence of px broadening of J/ip may not be a good signal 
for the deconfining phase transition. 

It is well known that in pA and AA collisions, the sec- 
ondary hadrons generally shows a pt broadening. The 
natural basis for this broadening is the initial state par- 
ton scatterings. For J/tp's, gluon fusion being the domi- 
nant mechanism for cc production, initial state scattering 
of the projectile/target gluons with the target/projectile 
nucleons causes the intrinsic momentum broadening of 
the gluons, which is reflected in the pt distribution of 
the resulting J/ip's. Parametrizing the intrinsic trans- 
verse momentum of a gluon, inside a nucleon as, 



f(q T ) ~ exp(-q T / < q% >) 



(1) 



momentum distribution of the resulting J/ip in NN col- 
lision is obtained by convoluting two such distributions, 



Inn (pt) ~ exp{-p T / < p\ >nn) 



(2) 

where < p\ >nn—< q\ > + < q\ >■ In nucleus- 
nucleus collisions at impact parameter b, if before fu- 
sion, a gluon undergo random walk and suffer N num- 
ber of subcollisions, its square momentum will increase 
to q T — ► Qt "b -^Oi $o being the average broadening in 
each subcollisions. Square momentum of J/ip then easily 
obtained as, 



< Pt >ab ( 5 ) =< Pt >nn +6 N AB (b) 



(3) 



where Nab (b) is the number of subcollisions suffered by 
the projectile and target gluons with the target and pro- 
jectile nucleons respectively. 

Average number of collisions A^s(b) can be obtained 
in a Glauber model. At impact parameter b, the posi- 
tions (s, z) and (b — s, z') specifies the formation point 
of cc in the two nuclei, with s in the transverse plane and 
z, z' along the beam axis. The number of collisions, prior 
to cc pair formation, can be written as, 



N(b,s,z,z') = (JgN dz A p A (s,z A ) (4) 

J —OQ 

+<T gN / dz B pB(b- s,z') 



1 



where a g N is the gluon-nucleon cross section. Above ex- 
pression should be averaged over all positions of cc for- 
mation with a weight given by the product of nuclear 
densities and survival probabilities S, 

//>oc />oo 
d 2 s / dzp A (s,z) dz'p B (b - s,z') x (5) 
J — oo J — oo 

S(b,s,z,z')N(b,s,z,z')/ / d 2 s / dz/^s, z) x 

J J —oo 

/oo 
dz'p B (b — s, z')S(b, s, z, z') 
-oo 

Finally, corresponding quantity at fixed transverse en- 
ergy Et is obtained as, 



Nab(Et) = / d 2 bP(b,E T )a AB N AB {b)/ 



d 2 bP(b, E T )a AB 



(6) 



/ 



where a AB is the inelastic cross section for AB collisions. 
P(b, Et) is the Ex — b correlation function. We have used 
the Gaussian form for the Et — b correlation, 

P(b, E T ) cx exp(-(E T - qN p (b)) 2 /2q 2 aN p (b)) (7) 

where N p (b) is the number of participant nucleons at 
impact parameter b. a and q are parameters related to 
dispersion and average transverse energy. For Pb+Pb 
collisions the parameters are, a=1.27 and g=0.274 GeV 
[3]. 

Survival probability 5(6, s, z, z') in Eq.5 is model de- 
pendent. We have calculated it using our QCD based 
'unconventional' nuclear absorption model [6,10]. Briefly, 
J/ip production is assumed to be a two step process, (a) 
formation of a cc pair, which is accurately calculable in 
QCD and (b) formation of a J/ip meson from the cc pair, 
which is conveniently parameterized. The J/ip cross sec- 
tion in AB collisions, at center of mass energy yfs was 
then written as, 



a J ^{s)=K 



a.b 



dq 2 



Cab- 



Q 2 



dx F <p a/A {x ai Q 2 ) (8) 



h/B(xb,Q 2 ) 



x F c5 ^j /i: (q 2 ), 



where ^ a b runs over all parton flavors, and Q 2 = q 2 + 
Am 2 . The K factor takes into account the higher order 
corrections. The incoming parton momentum fractions 
are fixed by kinematics and are x a = {\J x 2 F + AQ 2 /s + 
x F )/2 and x b = (y/ ' x 2 F + 4Q 2 /s - x F )/2. a ab ^ cS 
are the subprocess cross section and are given in [11]. 
Fcc->j/ip(q 2 ) is the transition probability that a cc pair 
with relative momentum square q 2 evolve into a physical 
J/ip meson. It is parameterized as, 

F ce ^j,4q 2 ) - N J/4 ,6(q 2 )6(4m' 2 - Am 2 c - q 2 ) (9) 



(1 



Am 



/2 



Am 2 



In a nucleon- nucleus /nucleus- nucleus collision, the pro- 
duced cc pairs interact with nuclear medium before they 
exit. It is argued that the interaction of a cc pair with 
nuclear environment increases the square of the relative 
momentum between the cc pair. As a result, some of 
the cc pairs can gain enough relative square momentum 
to cross the threshold to become an open charm meson. 
Consequently, the cross section for J/ ip production is re- 
duced in comparison with nucleon-nucleon cross section. 
If the J ftp meson travel a distance L, q 2 in the transition 



probability is replaced to q 2 



q 2 + e 2 L, 



e being the 



relative square momentum gain per unit length. Param- 
eters of the model were fixed from experimental data on 
total J/ip cross section in pA/AA collisions. It is thus 
essentially a parameter free calculation for Pb+Pb colli- 
sions. 

Fluctuations of Et at a fixed impact parameter plays 
an important role in J/ip suppression in Pb+Pb colli- 
sions. The 2nd drop in the J/tp over Drell-Yan ratio at 
100 GeV is due these fluctuations only. Fluctuations of 
Et at a fixed impact parameter also affect the average 
number of collisions N ab (Et)- As will be shown later, 
it plays an important role in explaining the NA50 pt 
broadening data. We have taken into account the Et 
fluctuations by the replacement, 



N AB {b) -> E T / <E t > (b)N AB (b). 



(10) 
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FIG. 1. (a) NA38 data on Et dependence of J /ip suppres- 
sion in S+U collisions. The solid line is a fit obtained in the 
'unconventional' nuclear absorption model, (b) NA38 data 
on pt broadening in S+U along with the fit obtained in the 
'unconventional' nuclear absorption model. 

Pt broadening of J/ip's in AA collisions depends on 
two parameters, (i) < p\ >nn, the mean squared trans- 
verse momentum in NN collisions, a measurable quantity 



2 



and (ii) the product of the gluon-nucleon cross section 
and the average parton momentum broadening per colli- 
sion, TgiySo- Since gluons arc not free, the second quan- 
tity is essentially non measurable. We obtain a g N8o from 
a fit to the NA38 pr broadening data [12] in S+U colli- 
sions at 200 GeV/c. <p\ >nn at corresponding energy 
is known from NA3 experiment, < p\ >nn— 1.23 ±0.05 
[13]. The Et — b correlation parameters, a and q for 
S+U collisions are, a=3.2 and g=0.74 GcV [14]. To 
show that the present 'unconventional' nuclear absorp- 
tion model also reproduces the centrality dependence of 
J/tp over Drell-Yan ratio in S+U collisions, in Fig. la, we 
have compared our results with the experimental data. 
We have neglected the effect of Et fluctuations in S+U 
collisions. The agreement between data and theory is 
good. In Fig. lb, NA38 experimental data on the Et de- 
pendence of pt broadening are shown. The solid line 
is a fit to the data, obtained with < p\ >nn= 1-23 
(fixed) and a g NS = 0.442 + 0.056. Value of a g ]ySo agrees 
closely with the value obtained by Kharzeev et al [9] in 
the conventional nuclear absorption model and also with 
the value obtained in the comover model [15]. 
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FIG. 2. Experimental < p T >nn as a function of cm en- 
ergy along with the fit with Eq.ll. 

< p T > nn increases weakly with energy. To obtain 

< Pt >nn f° r Pb+Pb collisions at 158 GeV/c, we have 
fitted the existing experimental data [13,16-18] with log- 
arithmic energy dependence, 

< Pt >nn= a + feln^/s (11) 

In Fig. 2, experimental data along with the fitted curve 
is shown. From the above parameterization, we obtain, 

< P T >nn =1-15 GeV 2 , for Pb+Pb collisions at CERN 
SPS. As we intend to predict pt broadening at RHIC en- 
ergy, < p T >nn at RHIC energy (-y/s=200 GeV) is also 
obtained from the above parameterization. At RHIC en- 
ergy, < p T >nn =2.45 GeV 2 . However, we must warn 
our reader to treat the above number with caution. The 



experimental data being limited to 60 GeV only, extrap- 
olation to RHIC energy is unreliable. 
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FIG. 3. (a) NA50 data on the centrality dependence of the 
ratio, J/tp over Drell-Yan, is compared with the 'unconven- 
tional' QCD based nuclear absorption model (solid line) and 
QGP based threshold model (dashed line), (b) NA50 data on 
p T broadening of J/ip's in Pb+Pb collisions. The solid and 
dashed lines are the prediction in the 'unconventional' nuclear 
absorption model, with and without the effect of Et fluctua- 
tions on NAB(b). The prediction in the QGP based threshold 
model [3], with and without the Et fluctuations are shown as 
the dash-dot and dash-dot-dot lines. 

For completeness purpose, in Fig. 3a, we have shown 
the NA50 data [2] on the centrality dependence of the 
ratio, J/tp over Drell-Yan. The solid line is the fit ob- 
tained to the data in the 'unconventional' QCD based 
nuclear absorption model. The data are well explained 
in the model. In Fig. 3a, we have also shown the J/tp 
suppression obtained in the QGP base threshold model 
[3]. In the threshold model, in addition to the conven- 
tional nuclear absorption, an anomalous suppression is 
included such that all the J ftp's are suppressed above 
a critical density n c . The dashed line is obtained in the 
threshold model with n c —3.7 fm 2 . [3]. It also gives satis- 
factory description to the data. Centrality dependence of 
J/tp suppression could not distinguish different natures 
of absorption. 

In Fig. 3b, we have compared the centrality dependence 
of pt broadening in the model with the NA50 experi- 
ment [8]. We have used < p T > NN —1.15 GeV 2 and 
Cgiv<>o=0.442. The solid and dashed lines are the pt 
broadening with and without the effect of Et fluctua- 
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tions on A/^s(6). When the effect of Et fluctuations 
is not taken into account, the pt broadening continues 
to increase with Et till 100 GeV (the knee of the Et 
distribution). Thereafter, pt broadening decreases. The 
behavior is unlike the pt broadening in the 'conventional' 
nuclear absorption model, rather more like the behavior 
in a QGP [9]. In the conventional nuclear absorption 
model, at large Et, Pt broadening saturates while in a 
QGP it decreases. Indeed, the different centrality depen- 
dence of pt broadening in nuclear and in QGP medium 
led Kharzeev et al [9] to suggest 'decreasing p T broaden- 
ing at large £t' as a signal of QGP. We find that 'un- 
conventional' nuclear absorption model also produces a 
'decreasing pt broadening at large Et'- At large Et, 
J ftp's are largely suppressed and < p T > decreases. De- 
crease of < p T > at large Et can not be considered as a 
signal of deconfincmcnt phase transition. 

Centrality dependence of < p\ > at large Et is 
changed when the effect of Et fluctuations, on average 
number of gluon-nucleon collisions, is taken into account. 
The decreasing trend of < p\ > beyond 100 GeV is 
changed to an increasing trend (the solid line). Effect of 
Et fluctuations essentially increases the average number 
of gluon-nucleon collisions at a fixed impact parameter 
and counter balance the large suppression effect beyond 
the knee of the Et distribution. NA50 data also shows an 
increasing trend beyond 100 GeV. The model reproduces 
the data very well (within 2%). 

We have also analysed the data in the QGP based 
threshold model [3]. Hucfncr et al [20] calculated pt 
broadening in the threshold model (Fig. 2 in ref [20]) 
and found that the NA50 data could not be fitted in 
the model. They did not take into account the effect of 
Et fluctuations on the average number of gluon-nucleon 
collisions, which we have seen, is important for explain- 
ing the data. In Fig. 3b, the dash-dot-dot line is the pt 
broadening in the threshold model, without the effect of 
Et fluctuations on Nab- As in [20] the model can not 
explain the data. However, at low Et, the it predict pt 
broadening in close agreement with the 'unconventional' 
QCD based nuclear absorption model. It is expected, as 
in peripheral collisions (low Et) QGP is not produced 
and J ftp's are suppressed in nuclear medium only. When 
the effect of Et fluctuations on NabQ)) is taken into ac- 
count (the dash-dot line), even though the difference be- 
tween the theory and experiment is lessened, the model 
still fails to explain the data. At large Et the model pre- 
dict 4% less pt broadening, also the upward trend beyond 
100 GeV is not reproduced. However, we note that the 
Pt broadening in the threshold model and in the nuclear 
absorption model are very close to each other, the differ- 
ence is less than within 3% or so. Since the NA50 data 
are very accurate, it can distinguish the small difference 
between the two models. We may mention that NA50 pt 
broadening data were analysed by Armesto et al [15], in 
the comover model. The comover model also gives very 
good fit to the data, only the increasing trend after Et= 
100 GeV, could not be reproduced. They also neglected 



the effect of Et fluctuations on the average number of 
gluon-nucleon collisions, which give that tendency. 

We now give prediction for pt broadening at RHIC 
energy. For RHIC energy, parameters a and q in the 
Et — b correlation arc taken 1.97 and q=0.46 GeV 

[7]. They were obtained by fitting rescaled Et distribu- 
tion data in Pb+Pb collisions at CERN SPS. At RHIC 
energy, the so called hard component, which is propor- 
tional to the number of binary collisions, appear. Model 
dependent calculations indicate that the hard component 
grows from 22% to 37% as the energy changes from 56 
GeV to 130 GeV [21]. In our calculation, we have used 
37% hard scattering component. In Fig. 4a, the predicted 
centrality dependence of the ratio of Jf tp over Drell-Yan 
ratio is shown. The QCD based nuclear absorption model 
and the QGP based threshold model, both predict nearly 
the same J/tp suppression. Different nature suppression 
is not distinguished. 
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FIG. 4. (a) Predicted centrality dependence of J/tp over 
Drell-Yan ratio at RHIC. The solid and dashed lines are the 
prediction in the QCD based nuclear absorption model and 
the QGP based threshold model, (b) Predicted centrality 
dependence of pr broadening of J/tp's. 

In Fig. 4b, the predicted pt broadening at RHIC are 
shown. We have used < p T >nn =2.45 (GeV 2 ) and 
Cgjv<5o=0-442. We again warn that these numbers are 
approximate only, The solid and the dashed lines are the 
predicted pt broadening obtained in the 'unconventional' 
nuclear absorption model and in the QGP based thresh- 
old model, both with the effect of Et fluctuations on 
Nab included. Both the models give similar centrality 
dependence for the pt broadening. After the initial rise 
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with Et, large suppression forces the px broadening to 
decrease, till 180 GeV( the knee of the Et distribution). 
Beyond the knee, large suppression is counter balanced 
by the effect of Et fluctuations on the average num- 
ber of gluon-nucleon collisions and pt broadening nearly 
saturates. The predictions in two models closely agree 
(within 1.5%) with each other. Given the uncertainty in 
< Pt >nn and a g N8o, such small difference may not be 
distinguished experimentally The results suggest that at 
RHIC energy, centrality dependence of pt broadening of 
J/ip, like the centrality dependence of J/ip, may not be 
able to distinguish a de-confinement phase transition. 

To summarize, we have analysed the NA50 data on pt 
broadening of J/ip's. The data were analysed in the QCD 
based, 'unconventional' nuclear absorption model [6] and 
in the QGP based threshold model [3]. It was shown that 
the Et fluctuations at a fixed impact parameter plays an 
important role in explaining the pr broadening data. If 
the effect is not incorporated, both the models predict 
a decreasing pt broadening at large Et, contrary to the 
experiment. When the effect is included, the 'unconven- 
tional' nuclear absorption model could explains the data 
without any free parameter. The QGP based threshold 
model could not explain the data. At large Et, it pro- 
duces 4% less pt broadening, also the increasing trend 
beyond 100 GeV is not reproduced. The analysis also 
indicate that, the 'visibly decreasing pt broadening at 
large Et ', can not be considered as a probe of the de- 
confinement phase transition. "Unconventional' nuclear 
absorption also produces a decreasing p T broadening at 
large Et- We have also obtained prediction for centrality 
dependence of pt broadening at RHIC energy. At RHIC, 
both the models predict very similar pt broadening. We 
conclude that pt broadening of J/^'s can not probe the 
deconfinement transition at RHIC energy also. 
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